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We study the flow close to an advancing contact line in the limit of small capillary 
number. To take into account wetting effects, both long and short-ranged contribu- 
tions to the disjoining pressure are taken into account. In front of the contact line, 
there is a microscopic film corresponding to a minimum of the interaction potential. 
We compute the parameters of the contact line solution relevant to the matching to a 
macroscopic problem, for example a spreading droplet. The result closely resembles 
previous results obtained with a slip model. 

I. INTRODUCTION 



Moving contact lines are encountered in a great number of flow problems, such as spread- 

n n 



ing of liquid drops jlj, dewetting of liquid films |2j], coating 3], and sloshing |4(. It was 
discovered by Huh and Scriven [5| that the viscous dissipation in the fluid wedge bordered 
by a solid and a fluid-gas interface is logarithmically infinite if the standard hydrodynamic 
equations and boundary conditions are used {(J. Thus continuum hydrodynamics does not 
describe the spreading of a drop on a table. Instead, some microscopic length scale must be 
introduced into the problem. 

As a model problem, let us consider the spreading of a viscous drop on a flat substrate. 
Typical spreading speeds are so small [1] that the bulk of the drop is almost unaffected by 
viscous shear forces. Hence the drop has the shape of a spherical cap, except in a small 
region around the contact line [t|. If one extrapolates this spherical cap solution to the 
contact line, it meets the solid at a well-defined angle, called the "apparent" contact angle 
9 ap . If for simplicity one assumes that the drop is thin, its radius R is related to 9 ap by 

e ap = AV/(7cR 3 ), (1) 

where V is the volume of the drop. 
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However, near the contact line the shear rate is of order U / h, where U is the contact line 
speed and h the local thickness of the fluid film. Near the contact line viscous forces become 
very large, and strongly bend the interface. A dimensionless measure of this viscous bending 
is the capillary number Ca = 7)U/j, representing a ratio of viscous to capillary forces, with 
77 the viscosity and 7 surface tension. As we will show below, within the approximation we 
adopt here, the slope h! of the interface as function of the distance x from the contact line 

n 

has the form |8|] 

h' 3 {x) = 9 3 e +9Caln{x/L), (2) 

where 9 e is the equilibrium contact angle and L a microscopic length scale. As illustrated 
in Fig. ^ we have adopted a coordinate system in which the contact line is at rest. The 
local description (J2J) applies for x/L 3> 1, i.e. at a distance from the contact line where 
microscopic details no longer matter. 

The distinguishing feature of (j2J) is that the curvature vanishes for x/L — > 00. This is 
a necessary condition for the local profile (J2J to be matchable to the spherical cap solution 
that makes up the bulk of the spreading drop |9( . The details of this matching procedure 
have been given in loj . the result being 

9% = 9 3 e + 9J&7/7 \n[R/{2e % L)\, (3) 

where e = 2.718281. . .. Together with (JSj) is evidently a differential equation for the 
radius of the spreading drop. For 9 ap ^> 9 e equations (JTJ), Q reproduce Tanner's spreading 
law Q R = At 1 / 10 , neglecting logarithmic corrections in time t. To find an explicit expression 
for A, it remains to know the length L. In this paper, we are going to compute L for a 
model that includes both long and short-ranged interactions in the interface potential Q]. 
This model has recently become popular for the numerical treatment of moving contact line 



problems 
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To find L, (J2J) has to be continued to the contact line, where microscopic effects come 
into pla y. P revious calculations have done that for the case of fluid slip over the solid 
surface |l3L Il4 L which relieves the contact line singularity. In the simplest case of a Navier 
slip condition described by a slip length A, the result is L = 3A/ (e9 e ). In we have 



extended this calculation to higher orders in the capillary number. However, corrections are 
found to be small in a regime where the underlying lubrication description is still expected 
to be valid ^(|. Apart from the slip length, an angle has to be specified at the contact line, 
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U 

FIG. 1: A cartoon of the contact line. In a frame of reference in which the contact line position is 
stationary, the solid moves to the right with velocity U. There is a very thin film of thickness hf 
in front of the contact line. 



which is often taken to be the equilibrium contact angle. This assumption implies that the 
total dissipation near the contact line is dominated by viscous effects, rather than dissipation 
localized at the contact line ~v\ . 

Here we adopt a model that builds in the equilibrium properties in a more rational way, 
w including the interface potential into the description. Both the equilibrium contact angle 
181 ] and the equilibrium film thickness h eg are determined by the interface potential. Within 
the model, even the "dry" substrate is covered by a thin film, corresponding to the minimum 
of the interface potential. The presence of this film thus formally eliminates the contact line 
singularity, h eq replacing the slip length A as the cut-off length. Of course, we do not claim 
that this is a true resolution of the contact line problem. The thickness h eq is often below 
the thickness of a single molecule, and even a monomolecular layer is not strictly describable 
by a continuum theory. 

Nevertheless, we believe that it is interesting to investigate the interplay between the 
interface potential and viscous forces. This has first been done by de Gennes, Hua, and 
Levinson Q], but only taking into account the long-ranged part of the potential. As a 
result, the equilibrium contact angle could only be worked in in an ad- hoc fashion, as one 
needs the full potential to define it. We will see below that our results are in line with the 
results obtained before [l^. The calculation in is based on a simple energy balance, 
rather than the systematic expansion performed here. The very recent work [21] treats both 
the advancing and the receding contact line in a manner very close to ours. 

Our paper is organized as follows. After introducing the model description, we recall the 
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case of a static contact line, relating the equilibrium contact angle to the interface potential. 
We then outline how the parameter L of (J2J) may be found in an expansion in the capillary 
number Assuming a particular form of the interface potential, we then solve the first 
order problem explicitly. Finally, we compare to other forms of the interface potential as 
well as to previous work. 



II. LUBRICATION DESCRIPTION 



For simplicity, we perform our calculations within the framework of lubrication theory, 
;hus limiting ourselves to the case of small contact angles, as well as small capillary number 



_22| . Experiment shows that this approximation performs reasonably well up to a capillary 
number of 0.1 Q]. The lubrication equation reads ^| 

3rjh t = -\h 3 ^h xx + U(h)) x ] , (4) 



where h(x,t) is the thickness of the fluid film and n(/i) is the disjoining pressure [18(. The 
origin of (jlj) is a viscous shear flow, driven by the gradient of the pressure p = —^h xx — U(h). 
The first term is the usual Laplace pressure, proportional to the curvature of the interface, 
while the disjoining pressure H(h) is given by H(h) = dV/dh, where V(h) is the effective 
interface potential of a flat film of thickness h [2] . Thus as soon as h is larger than the range 
of all the interactions between particles, H(h) can safely be neglected. However, when h is 
of the order of a few nanometers, the disjoining pressure becomes relevant. 

To describe an advancing contact line (cf. Fig. QJ, it is convenient to pass into a frame 
of reference that moves with the contact line speed U: 

h(x,t) = h(x + Ut), (5) 

giving 

3Cah x = - [h\h xx + U(h)/j) x ] x . (6) 
Integrating once one finds that 

-[h xx + U(h)/j] x , (7) 



2>Ca{h - h f ) 



h 3 

where hf is the (yet unknown) film thickness ahead of the moving contact line. 
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III. STATICS 

It is instructive to look first at the well-known static case Ca = 0. Integrating (J7J) once 
more one obtains 

Po = ~h xx - n(/i)/ 7 , (8) 

where Po is the (constant) pressure in the film (neglecting gravity). We are considering 
a situation where the film is in contact with a large reservoir (for example a drop) with 
negligible pressure, hence Po = 0. Thus in the film we must have U(h eq ) = (corresponding 
to a minimum of the interface potential), which defines the equilibrium film thickness h eq . 
Now (jBJ) can easily be solved by putting g(h) = h x (x), giving 

^ = -2U(h)/ r (9) 

Integrating (jHJ), we obtain the standard expression [lsl ] 

e 2 e = -2 n(c)/ 7 dC (io) 

for the equilibrium contact angle, which in the lubrication approximation is to be identified 
with the slope of the interface: 9 e = tan(/i x (oo)) fa h x (oo). By integrating to infinity, we 
imply that the macroscopic scale on which 9 e is defined is much larger than h eq . 

To be more specific, the disjoining pressure has a long-ranged attractive and a short- 
ranged repulsive part: 

A R 

n(/i) = — (n) 

The repulsive interaction keeps the film thickness from collapsing to zero. The form of the 
attractive part is rather universal pj, A being known as the Hamaker constant. The most 
popular choice for the repulsive part is a power law with a = 9, which is motivated by 
the form of the Lennard- Jones interaction. Recently, enormous progress has been made in 
determining the constants in (fTT|) for some systems 2]. However, the experiments are not 
sufficiently accurate to determine the value of the exponent a j^. For some of the explicit 
results to be reported below we are going to choose another value, a = 5, to be able to 
perform our calculations analytically. Using the specific form of (jlljl . one easily finds that 

heq = (B/A) in "~ 3 \ 01^-1^. (12) 
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FIG. 2: Numerical solutions for the rescaled interface slope -ff (£)• The full line is the static solution 
()16|) . the dashed line a solution of (|17|) for 5 = 0.1 and a = 5. To expand the horizontal range, 
a logarithmic scale has been chosen, with £ = shifted somewhat to the left of the contact line 
position. 

To compute the profile, it is useful to introduce new variables, which are scaled to the 
equilibrium thickness h eq of the film: 



h(x) = h eq H(£), i = xdJK 



eg • 



(13) 



Equation Q then becomes 



H 



a — 1 , . 
2 — 



a-3\H 3 H a , 

To make further progress, we specialize to a = 5, in which case we simply have: 

H 2 -l 



(14) 



H' 



H 2 



(15) 



This can be integrated to give the static interface shape 

«^ + 5 ln (fTl)- (16) 

where the left hand side can of course be shifted by an arbitrary amount. The slope of the 
static interface is shown in Fig|21 To the right of the contact line the slope asymptotes to 
1, corresponding to the equilibrium contact angle. 
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IV. PERTURBATION EXPANSION 

Now we turn to the problem of a moving contact line. In the scaled description (|13jh (J7J) 

becomes 

S(H - H f ) 



H 3 



_ H „ + a-lfl 1 



a -3 \H 3 H a J\ ' ^ 
where 5 = ?>Ca/9 3 is the rescaled capillary number. In the limit of small h eq the boundaries 
of the system are pushed out to £ = ±oo, and the boundary conditions become 

H(-oo) = H f , H'(-oo) = 0, H"(oo) = 0. (18) 

The first two conditions correspond to the assumption that the liquid forms a film of constant 
thickness ahead of the contact line. We will see below that it deviates slightly from the 
equilibrium thickness if the contact line is moving. The third boundary condition says that 
the curvature far away form the contact line is vanishingly small compared to the typical 
curvature near the contact line, which is l/h eq Q|. 

We are going to solve (fTT|) in a perturbation expansion in S, following a procedure adopted 
before Of particular interest is the behavior of the solution for large £, which corre- 
sponds to (0). Namely, for H ^> 1 (fTTj) assumes the universal form 5/H 2 = H'", which has 
the asymptotic solution |24j 

tf'(0 = [351n(^o)] 1/3 , £»1. (19) 

This solution has vanishing curvature at infinity (as required by (fTBj) ). and only contains a 
single free parameter £ > to be determined by matching to the contact line. By comparing 
(HH and (J2J), one finds 

^ = ioe^\ (20) 

''■eg 

On the other hand, the full solution H(£) possesses a perturbation expansion in 5 around 
the static profile Hq(£) '■ 

H(0 = H (0 + 6H l (0 + O(6 2 ). (21) 

For large £, we have -£^(0 ~ 1, corresponding to the equilibrium contact angle. By com- 
paring this to (fTTH) . we find that ln(£o) has the following expansion: 

-31n(£ ) = ^ + cl + O(5). (22) 
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Substituting into (fT§j). we find that for large £ 

tf((0 = ln(0 + ci/3. (23) 

To compute L, we thus take the following steps: First, we solve the full problem (JT7jl 
perturbatively to obtain Hi(£). Then, analyzing Hi for large £, we obtain ci, which gives 
£o by virtue of ()22jl . Combining this with ()2()jl . we finally have 

L = ^ie" C1 / 3 . (24) 

V. EXPLICIT SOLUTION 

To first order in S, (|17|) becomes 



5 ^« = -^+ — (^+-^)+c' I (25) 

-oo ff 3 S 1 a - 3 V #o #o ' 



where we have integrated once, resulting in a constant of integration C. From now on we 
consider the special case a = 5, for which we can make use of the static solution i?o(0 given 
by (HSJ). 

The integral on the left-hand-side of (|23|) can be performed by exchanging the role of 
dependent and independent variables using (|15j) : 

« #o - 1 ^ / ,fr ° Irr , / 2^o 



, <%= ; zdH Q = In — . (26) 

Hi ^ Ji H (H + 1) \Ho + l) V ' 

The limit of (|26|) for large £ is ln(2), hence taking the same limit in (J25j) yields C = ln(2) 
for the constant of integration. Now considering the opposite limit of £ — > — oo, and using 
Fo(-oo) = 1, one finds #i(-oo) = -ln(2)/4. 

To solve (|23j). it is useful to rewrite the entire equation using H as the independent 
variable. To avoid cumbersome expressions, we denote Ho by the symbol (. Thus ()25)1 turns 
into: 

F K> - ln (cTl) ^ -™« + 2 ™< (? " ' < 27 > 

Remarkably, this equation can be solved exactly by noticing that two fundamental solutions 
are 



KD 16C 5 -5QC J + 30C C J -i ln fC-i\ and H m_C_lz 
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FIG. 3: A comparison with simulation. The full line is the leading-order result of our calculation 
(jSU, the dashed line is the numerical result, valid to all orders in 5. 



which we found using Maple. Thus a general solution of ()27j) is 



b x - ( H^F/WdQ' +H[ 2) b 2 + [ C H[ 1] F/WdC 

J2 J2 



(29) 



Hi = HI 1 

where W is the Wronskian. 

The limit ( — > 1 corresponds to the thin film. From the condition that Hi has to remain 
finite in this limit, one finds 

h = - J* H^F/WdC = 3 ln(3)/16 - ln(2)/4, (30) 

since Hi — > oo for ( — > 1. As shown in the Appendix, the other constant of integration 
&2 is determined by the terms of order (° as ( — > 1. In the limit of ( — > oo, on the other 
hand, one is approaching the bulk fluid, for which we find f/j ~ 16^ and Hi ~ 1, so a 
straightforward analysis of (|29|) yields 

Hi(C) = C(ln(C) - 21n(2)) + O(m(0). (31) 



VI. RESULTS AND DISCUSSION 

Now we are in a position to calculate the constant Ci appearing in 
have C = H w ^ for large £, and thus 



. From ()16|) we 



^(0 = (ln(0-21n(2) + l) 



(32) 
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in this limit. We conclude that Hi indeed has the asymptotic form (J2~3"j) we anticipated, and 
we can identify 

ci = 3-61n(2). (33) 



Using (j24j) . we now have 



Ah 

^ (34) 



which is the central result of this paper. 

The result ()34j) can of course be tested by comparing with a numerical solution of the full 
equation (|17j). A linear analysis around the film thickness H = Hf reveals an exponentially 
growing solution 

#(0 = #/ + eexp(70, (35) 

where 7 = 2 + 0(5). Any small perturbation of the constant solution H = Hf, H' = 0, H" = 

_growth of the form ()35|). As £ — > 00, the solution generically 



will thus lead to an initial 
tends to a finite curvature 



24j . Thus Hf has to be adjusted to find the unique solution 



which obeys the boundary condition ([18)1 at infinity. The asymptotics of this solution of 
course has to conform with (|19j) . 

However, the approach to this solution is very slow, as revealed by the full asymptotic 
expansion |25J] 

{°° h 1 

1+ 5(w|- (36) 

To be consistent with (j!9|) . the coefficient hi was chosen to vanish, since it would lead to 
a redefinition of £0 To obtain £ numerically, we fitted the numerical solution of (|17|) to 
(I36|) . using the first five terms of the expansion. In Fig. 01 we plot the numerical result for 
L over a wide range of 5-values. For reasonably small <5's, applicable to most experimental 
situations, the result is very well approximated by the present first order calculation. 

Our analytical approach has of course been limited to the case a = 5, which is non- 
standard. Using the numerical procedure described above, it is a simple matter to obtain L 
for arbitrary a. Fig0] shows the result of this calculation in the limit of small 5. As to be 
expected, the variation with a is not very strong. Large values of a correspond to a very 
hard core. 

Finally, it remains to compare our results to who only took the long-ranged part of 
the disjoining pressure into account. At the contact line, it was assumed that the solution 
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FIG. 4: The characteristic length L as function of the exponent a characterizing the potential. 
For a = 5 the numerical result agrees with (|34|). 



matches to the equilibrium contact angle. The result was reported in the form Lghl — 
a I (291), where 



is a length scale characterizing the range of van-der-Waals forces. Thus, using (112)1 the result 
of UM can be converted to 



which is essentially the same result as (jSJ), but with a different prefactor. In conclusion, for 
both a slip and the present thin film model, L is set by the respective microscopic length. 



I am grateful to Len Pismen for his input, and to the participants of the Thin Film 
workshop in Udine in 2005, organized by Serafim Kalliadasis and Uwe Thiele, for advice. 

APPENDIX 

Here we describe how to determine the remaining constant of integration 6 2 m (|29|). by 
comparing to the asymptotics (jH5]) of the full solution as £ — ► — oo. Namely, as we have 




(37) 




(38) 
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shown above, 

H f = l-(\n(2)/A)5 + 0(5 2 ), (A.l) 
and it is straightforward to see that the exponent is 

7 = 2 + 7l 5 + 0(5 2 ), 7 i = 91n(2)/4-l/8. (A.2) 

Thus at zeroth order in S one finds ( = 1 + eexp(2£). On the other hand, the full static 
profile (|TH|) gives 2(£ — 1 + ln(2)/2) = ln(£ — 1) + 0(( — 1). Thus by comparing the two 
profiles one identifies ln(e) = ln(2) — 2. 

Expanding (J3l)j) to next order in 5 leads to 

E x = - ln(2)/4 + 7l (f - 1 + ln(2)/2) exp[2£ - 2 + ln(2)] = (A.3) 
- ln(2)/4 + 7l ln(C - 1)(C - 1) + 0(C - l) 2 - 

Thus in the limit of £ — ► 1, (J2HJ) must have the same form as (|A.3j) . The integrals in (J2HJ) 
can be performed using Maple, and in the limit they give 

H x = - ln(2)/4 + (a + 7l ln(C - 1)) (C - 1), 

which matches (|A.3)1 if a = 0. From this requirement we finally get 

b * = " 1? ln ( 3 ) " 3 ( 2 / 3 ) + 3/2 di/pff (4/3) + (A.4) 

it; 7c ioc 

— dilog (3) + — (ln (3)) 2 + — ln (2) - 
16 y v ; 32 v v 48 V ; 

97 1 ^ 1 3 

^ ln (2) ln (3) - ^ (ln (2)) 2 + ^ vr 2 = 0.359777 .... 
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